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Abstract 



^ Grand canonical ensemble of small vortex loops emerging in the London limit of the effective 

fSl ' Abelian-projected theory of the S'?7(3)-gluodynamics is investigated in the dilute gas approxima- 

^ \ tion. An essential difference of this system from the SU (2)-case is the presence of two interacting 

I gases of vortex loops. Two alternative representations for the partition function of such a grand 

O ■ canonical ensemble are derived, and one of them, which is a representation in terms of the integrals 

I over vortex loops, is employed for the evaluation of the correlators of both kinds of loops in the 

0^ ! low-energy limit. 
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In a recent paper [Q, the grand canonical ensemble of small vortex loops, existing in the 
Abelian Higgs model, have been investigated. Such loops are nothing else but the 4D analogue 
of the vortex dipoles, which are present in the usual Ginzburg-Landau theory. In particular, it 
has been demonstrated that the summation over the grand canonical ensemble of small vortex 
loops leads to an effective Sine-Gordon type theory of the massive Kalb-Ramond field [@ (c/. 2D- 
and 3D cases studied in Refs. 0, 0). This field describes a dual vector boson, which therefore 
acquires an additional mass due to the Debye screening in the gas of vortex loops. Furthermore, a 
representation of the partition function of such a gas directly in terms of the integral over vortex 
loops as well as the related effective potential of those have been discussed. Such a representation 
then turned out to be useful for the evaluation of the bilocal correlator of vortex loops in the low- 
energy limit. These calculations demonstrate the importance of treating the topological defects 
in the Abelian Higgs model and Ginzburg-Landau theory as forming the ensembles, rather than 
individual ones. 

An interest to the study of vortex loops in the Abelian Higgs model is motivated by the fact 
that the dual Abelian Higgs model is discussed to be relevant to the description of confinement 
in the S't/(2)-gluodynamics 0. This agreement is based on the method of Abelian projections 
(for a recent progress see 0, |, ^, [T^, |TT], |T2|, [T^, 0, |T3|, for a review see |]TB|) and the so-called 
Abelian dominance hypothesis |T^, according to which off-diagonal (in the sense of the Cartan 
decomposition) degrees of freedom are inessential for confinement and can be disregarded. In 
the spirit of this hypothesis, the following partition function describing an effective [[/(l)]^ gauge 
invariant Abelian-projected theory of the realistic SU (3)-gluodynamics has been proposed [l^ Q 



Z = JvB^VxaVxlS (^E^a) 



X exp <— d X 



1 
4 



r -►2 



a=l 



Here, F^y = d^B^ — d^B^ stands for the field strength tensor of the magnetic vector potential 



Ida 



1,2,3, 



= S«) dual to the electric one = (A^, A^). Next in Eq. (§, Xa = |Xa| e 
are effective Higgs fields describing condensed magnetic monopoles, whose magnetic charge Qm is 
expressed via the QCD coupling constant Qqcd as gm = '^'^/gqcD- Finally in Eq. (|l|). 



£l = (l,0), €2 




£3 




stand for the so-called root vectors, which play the role of the structural constants in the algebra 
H,E±a = ^£aE±a- Here, the operators H = {Hi,H2) = (T^jTs) generate the Cartan subal- 
gebra, where from now on Tj = z = 1, . . . , 8, are just the S'[/(3)-generators. We have also 
introduced the so-called step operators E±aS (else called raising operators for positive a's and 
lowering operators otherwise) by redefining the rest (non-diagonal) S'f/(3)-generators as follows 



±1 



±2 



±3 



^ Throughout the present Letter, aU the investigations will be performed in the Euclidean space-time. 
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Clearly, these operators are non-Hermitean in the sense that (Ea) = E^a- For bookkeeping 
purposes, it is worth listing the remaining commutation relations, completing the Lie algebra, 
which read 



[E±a, E±b\ — T-^^abcEz^c 



— * 

and [Ea, = SabSaH. 



Notice also that due to the fact that the original SU (3) group is special, the phases of the three 

3 

magnetic Higgs fields are not independent and should obey the constraint 6*0 = 0. The latter 

a=l 

one has been imposed by the introduction of the corresponding 5-function into the functional 
integral on the R.H.S. of Eq. (|I]). 

Before proceeding with the study of the model (P, it is worth mentioning its certain feature 
owing to which this model is not quite adequate to the description of confinement in the real 
S'f/(3)-gluodynamics. Its essence is that the model (|I]) describes only the sector of the full yet un- 
known Abelian-projected theory of the SU (3)-gluodynamics, where antimonopoles are completely 
absent. Clearly, the expected full theory should contain the antimonopole sector as well. Possible 
interference between these two sectors is up to now unclear and should be clarified by further 
investigations. 

In what follows, we shall be interested in the study of the model ([T|) in the London limit, i.e., 
the limit of infinitely large Higgs coupling constant A. Analogously to the 5'f/(2)-case, in this 
limit the model under study allows for an exact reformulation in terms of the integral over closed 
Abrikosov- Nielsen- Olesen type strings [^] 0. In this limit, the radial parts of the Higgs fields can 
be integrated out, and the partition function takes the form 




exp 



d X 



1 



a=l 



(2) 



Next, the total phases 's of magnetic Higgs fields should be decomposed into the singular and 
regular parts, 9 a = + 0> S 11 (c/ also Refs. P, |l5l for the Sf/(2)-case). Here, 



^^™^'s describe a certain configuration of electric strings and are unambiguously related to their 
world-sheets S^'s according to the equation (see the above cited Refs.) 



2vrS»,(x) 



27r / da 



(3) 



This equation is just the covariant formulation of the 4D analogue of the Stokes theorem for 



the gradient of the field 9a, written in the local form. In Eq. 



(0 



— ™(a) 



(,^) is a vector 



parametrizing the world-sheet Ha with ^ = {^^, G [0, 1] x [0, 1] standing for the two-dimensional 
coordinate. 

Confining and topological properties of the model (0) have been studied in Refs. ||I2|, |13|, |n . 
This has been done by making use of the so-called path-integral duality transformation, elaborated 
in Refs. ^ for the usual Abelian Higgs model, which casts the partition function into the 
following form, 

Z = J (0^ VAlVhl^ e^v\-j d\ ^ (i/^^,,) ' + 

19], it is this limit ol the Abelian-projected theories, in which they 



^Notice that according to the lattice data 
reveal properties similar to the real QCD. 
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with Af, 



EaAu. Here, 



+ 9n 



Tja 



(4) 



+ duh^^ stands for the field strength tensor of the 



antisymmetric tensor field /i^^, (the so-called Kalb-Ramond field 0). The integration over this field 
came about via some constraints resulting from the integration over O^^^^s, whereas the integration 
over ^^™s'g iig^g transformed into the integration over x|^°'*(^)'s by virtue of Eq. (0) (Notice that since 
in what follows we shall be interested in effective actions rather than the integration measures, the 
Jacobian appearing during the change of the integration variables, 6'^™^ ^l^HO^ which has been 
evaluated in Ref. will not be discussed below and is assumed to be included into the measure 
Vx^i^).). Also, due to this one-to-one correspondence between 6'^™s's and S^'s, the constraint 
imposed by the 5-function on the R.H.S. of Eqs. (|l|) and (^ has gone over into the constraint 
imposed by the 5-function on the R.H.S. of Eq. (^, which relates the world-sheets of three types 
to each other, making only two of them really independent. 

The aim of the present Letter is to treat Abrikosov-Nielsen-Olesen type strings in the model (^ 
in the sense of the grand canonical ensemble of small vortex loops, rather than as individual {i.e., 
noninteracting) ones. To understand why one might expect in this case the appearance of some 
nontrivialities w.r.t. the Abelian-projected S'f/(2)-gluodynamics, let us begin with considering 
noninteracting vortex loops. This can be done by gauging the field away from Eq. (^) by 

performing the hypergauge transformation /i^^^ — > /i^^ — {ptJ.^u ~ "^j^^/I) ^'^'^ subsequent 
integration over the Kalb-Ramond fields (see the first paper from Ref. [|15| for details of this 
integration). The result has the form 



y\ 



where m = y/SgmV is the mass of the fields and B^, which they acquire due to the Higgs 
mechanism, and Ki stands for the modified Bessel function. Finally, one of the three world- 
sheets, for concreteness a;|f''(^), can be integrated out, which yields 



Z = jVx^^\i)Vxf{Ox 



X ( 



exp -g^T^'V^ ( d'xd'y [s;,(x)Si,(y) + Si,(x)Sj,(i/) + ^l{x)^liy) 



Kilmlx — y\ 
\x - y\ 



■ (5) 



In order to proceed from the individual strings to the grand canonical ensemble of interacting 
vortex loops, one should replace E^j^(x), where from now on a = 1, 2, in Eq. (|^) by 

N 

k=i 

Here, n^^^^s stand for winding numbers, which we shall set to be equal ±1 (c/. Ref. [T[]) Perform- 
ing such a replacement, one can see the crucial difference of the grand canonical ensemble of small 



This is just the essence of the dipole approximation. 
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vortex loops in the model under study from that in the Abelian-projected SU (2)-gluodynamics [|T|. 
Namely, the system has now the form of two interacting gases consisting of the vortex loops of 
two kinds, while in the SU (2)-case the gas was built out of vortex loops of the only one kind. 

Analogously to that case, we shall treat such a grand canonical ensemble of vortex loops in 
the dilute gas approximation. According to it, characteristic sizes of loops are much smaller 
than characteristic distances between them, which in particular means that the vortex loops are 
short living objects. Then the summation over this grand canonical ensemble can be most easily 
performed by inserting the unity 



yiagas 



(6) 



into the R.H.S. of Eq. @ (with S^^ replaced by S^^**^) and representing the 5-functions as the 
integrals over Lagrange multipliers. Then, the contribution of vortex loops of each kind to the 
full grand canonical ensemble takes the following form 



Z 



v-iagas 



X exp 



g^r^v'Vs [ d'xdS [sl^ix)Sl^iy) + Sl,ix)Sl^iy) + Sl,{x)SlM 



Ki{m\x — y\) 
\x-y\ 



gas 



(7) 



After that, the desired summation is straightforward, since it technically parallels the one of 
Abelian-projected S'?7(2)-gluodynamics described in Ref. We have 



oo /-AT / N 



C 



.(1) 



X 



J2 exp 



TV 

L k=l 



r 



X 



X 



{the same term with the replacements (1) (2) and A^^ — * -^^z^} 



exp I 2C J d^y 



cos 



cos 



KM 



A2 



Here, the world-sheet coordinate of the /c-th vortex loop of the a-th type Q x^lf'{^) has been 
decomposed as x^^\^) = y^"^ + zl^\^), where the vector y^^^ = J d'^^xl^\^) describes the position 
of the vortex loop, whereas the vector zl^\^) describes its shape. Next, on the L.H.S. of Eq. (|), 
stands for a certain rotation- and translation invariant measure of integration over the 



shapes of the world-sheets of the vortex loops, and C oc e "is the so-called fugacity (Boltzmann 



*For brevity, we omit the Lorentz index. 
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factor of a single vortex loop of dimension (mass)^ with 5*0 denoting the action of a single loop. 
In Eq. (1^), we have also introduced the UV momentum cutoff A = J a' 



where a is a 

typical size of the vortex loop, and L is a typical distance between loops, so that in the dilute gas 
approximation under study a <^ L. Finally in Eq. I 



we have denoted 



M ■ The 



reader is referred to Ref. |1|] for details of a derivation of Eq. 

Note that the value of 5*0 is approximately equal to aa^, where we have estimated the area of 
a vortex loop as a^, and a stands for an analogue of the string tension for the loop, i.e., its energy 
per unit area. This energy can be evaluated from Eq. (|^) by virtue of the results of Ref. and 
has the form 



a 



2r]' / dH 



,2 Km) 



27rry^ln(A' 



(9) 



Here, we have in the standard way [T^ set for a characteristic small dimensionless quantity in the 
model under study the value which is of the order of the ratio of m to the masses of magnetic 

Higgs fields. Moreover, it has been assumed that not only — ^ 1, but also In — ^ 1, i.e., the 
last equality on the R.H.S. of Eq. is valid with the logarithmic accuracy. 

Next, it is possible to integrate out the Lagrange multipliers by solving the saddle-point equa- 
tions following from Eqs. (0) and (H), 



sm 




iK^ 
' 2C 



After that, we arrive at the following representation for the partition function of the grand canon- 
ical ensemble 



2grand = / 'T^S'^^v ^^P ' 



g^^^^?> I d'xd'y [sl^{x)Sl{y) + Sl^{x)SlM + Sl^{x)SlM 



X 



^ Ki{m\x-y\) ^ ^ 
\x-y\ 



r.1 



V 



s\ 



(10) 



which owing to Eq. is natural to be referred to as the representation in terms of the vortex 
loops. In Eq. ([10|), the effective potential of vortex loops reads 



V 



na 



E 



d^x I A' 



-2C 



In 



A2 
2C 



+ 



\ 



1 + 



'A2 



Ca 



'A2 
2C 



Ca 



+ 27rm 



It is further instructive to illustrate the difference of such a partition function of two interacting 
gases of vortex loops from the case of Abelian-projected SU (2)-gluodynamics by studying a related 

^It is natural to assume that the vortex loops of different kinds have the same fugacity, since different 6'^'"S's 
enter the initial partition function (0) in the same way. 
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representation in terms of a certain effective Sine-Gordon theory. This can be done by introducing 
the new integration variables iS^^, = ^ (^Sj^^, + Sf^^^ and iS^^, = | (^Sj^^, — Sf^Jj , which diagonahze 
the quadratic form in square brackets on the R.H.S. of Eq. Then Eqs. and (§) yield 

-z I d'xh;,S;,y (12) 

where we have denoted /i^^ = (^A^^, + A^^^j and /i^^, = A^j, — A^^,. The partition function of 
the desired Sine-Gordon theory can be obtained from Eq. (|T2p by making use of the following 
equality 



d X 



cos 



A2 




cos 



A2 



gmV'V3 f d^xd^ySt 



X 



Ki{m\x — y\ 



\x - y\ 



exp 



47]-2 



d X 



+ 



(c/. the R.H.S. with the quadratic part of the action of the Kalb-Ramond field on the R.H.S. of 
Eq. @ with the field gauged away). Substituting this equality into Eq. (|12D and performing 

the rescaling ^ /ij^^^, we arrive at the following representation for the partition function of the 
grand canonical ensemble of vortex loops in terms of the local Sine-Gordon theory, equivalent to 
the nonlocal theory (|10|), 



"'fj.u 



-2C 



cos 



TT 
A2 



+ cos 



7r 

A2 



V3h 



fJ,U 



h 



(13) 



As we now see, an essential property of the obtained Sine-Gordon theory, which distinguishes it 
from an analogous theory describing the grand canonical ensemble of vortex loops in the Abelian- 
projected S'f/(2)-gluodynamics [|I|, is the presence of two interacting Kalb-Ramond fields, while 
in the 5'?7(2)-case there was only one self-interacting field. Notice that upon the expansion of the 
cosines on the R.H.S. of Eq. (|l3l), it is straightforward to see that only the interaction terms of the 

/ . \ 2n / \2k 

type [h^j^y] [hjiiy] survive. In another words, despite of the mixing of the Kalb-Ramond fields 
in the arguments of the cosines, no terms linear in any of these fields appear in the action. In 
particular, the masses of both Kalb-Ramond fields. Mi and M2, can be read off from Eq. (|13|) by 
expanding the cosines up to the quadratic terms. The result reads = w? -t- = Q^rf, where 



mi = ^v^, m2 = 
magnetic charges Qi 



^y2C are the Debye masses, and we have introduced the corresponding 



Eqs. (|r^) and (|TTp can now be used for the evaluation of correlators of vortex loops, which due 
to Eq. (0), are nothing else but the correlators of S°j^Js. Those are calculable in the low-energy 



^This equality can straightforwardly be proved by mentioning that 9^^^^ = 0, which follows from the equation 

according to which d^S'^^^ 



of motion corresponding to its L.H.S. and Eq. 



0. 
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limit, S^Jj, <^ (, by considering the real branch of the potential (|n|), i.e., extracting from the 
whole sum the term with n = 0. The latter one has a simple parabolic form, and in the vicinity 
of its minimum (corresponding to the low-energy limit) the generating functional for correlators 



of S^j^^s reads 



Z 



ja 



X 



Ki{m\x-y\) 



\x - y\ 



A' 



+ / d'^x 



j+ 



where J^^ is a source of 5*^^,, and = J^^ ± J^^,. Such two Gaussian integrals can be calculated 
by virtue of the following equality 



VSa^ exp 



gmV^V^ I d'^xd'^ySf.^ix) ^^^^^^ ^ ^^^ S^^{y) + ^ / d'^^'^lu + / d'^xJ^^S^^ 



\x - y\ 



= exp 

and the result reads 



M2C 
'87r2A4 



d'xd'yJ,^{x)J,^{y) (dl - m') ^'^^'^"^ , ^'^ 



x-y\ 



Z 



ja 



-^grand GXp 



c 



87r2A4 



d'^xd'^y 



\x - y\ 



.2 2 



K,{M2\x-y\) 



\x - y\ 



The correlators of vortex loops following from this expression have the form 

c 



2A4 



and 



(ml 
2A4 



M2K,{M2 




47r2 


x\ 





25ix)-J2ml 

a=l 



-3 





N) 


47]-2 


x\ 





Ml Ki{Mi 




47r2 


x\ 





At this point, it is worth recalling that the original theory (|I]) is an effective theory at large 
distances [0, where the asymptotic behaviours of the obtained correlators read 



SUx)Sl{0)) = (Sl{x)Sl{0) 



_2 Malxl 
3 ^ 



X 2 
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This result illustrates how the vortex loops in the grand canonical ensemble under study are 
correlated to each other. Namely, their correlators decrease according to the Yukawa type law 
with the screening provided by the lightest of the two full masses, M2. 

In conclusion, we have demonstrated that the grand canonical ensemble of vortex loops in the 
effective Abelian-projected theory of the S'?7(3)-gluodynamics (being treated in the dilute gas ap- 
proximation) exhibits an essential property distinguishing it from the one of the Abelian-projected 
SU (2)-gluodynamics. Namely, it consists of two interacting subsystems, corresponding to two in- 
dependent types of strings, which emerge after the Abelian projection. An average over the shapes 
of the vortex loops with the most general rotation- and translation invariant integration measure 
leads to two alternative field-theoretical representations of such a grand canonical ensemble. First 
of them is a representation in terms of an effective Sine-Gordon theory of two interacting Kalb- 
Ramond fields (|T^) . It yields the (positive) contribution to the masses M^s of the Kalb-Ramond 
fields coming about from the Debye screening as well as to the magnetic charges QaS of these 
fields. The other representation given by Eqs. ([T0|) and (|ll]), which is the one in terms of the inte- 
gral over the vortex loops, is useful for the evaluation of their correlators. While such a calculation 
is difficult to perform exactly due to the complicated form of the effective potential of the vortex 
loops ([TT|), it turns out to be possible to perform it in the low-energy limit within an additional 
approximation when only the real branch of the potential is taken into account. As a result, the 
correlators of vortex loops have a Yukawa type asymptotic behaviours at large distances with the 
screening governed by the lightest of the two full masses of the Kalb-Ramond fields. 

It now looks reasonable to find field-theoretical representations for the grand canonical ensem- 
bles of vortex loops emerging after the Abelian projection in the general SU{N), N > 2, case. 
In that case, there appear _ 1 independent strings. Indeed, this is just the number of 

possibilities for the eigenvalues of a certain (adjointly transformed) operator, to be diagonalized 
during the Abelian projection, to coincide, minus one constraint imposed by the 5-function on 
the R.H.S. of Eq. (|1|) (with the sum over a going from 1 to A^). It will therefore emerge just this 
amount of interacting gases of vortex loops, leading to different Debye masses. The study of this 
system as well as its large- limit will be the topic of a separate publication. 
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